Transmission networks inferred from HIV sequence data

HIV in the UK in the 1980s was concentrated within men who have sex with men (MSM) and people who inject drugs (PWID) but heterosexual sex is now the most frequently reported risk behaviour. As these risk groups are associated with different virus populations, this is reflected in the subtype diversification of the UK epidemic, which was historically dominated by subtype B. I have made use of a national database of HIV sequences collected during routine clinical care, which also contains data on age, sex, route of exposure & ethnicity. The 2014 release of the UK HIV Drug Resistance Database contained data from over 60,000 patients. In this thesis, I first describe the development of novel tools that rapidly and automatically identify HIV clusters within phylogenetic trees containing tens of thousands of sequences because they represent transmission chains within the larger infected population. I use these tools to compare the HIV subtype B epidemics in the UK and Switzerland, which had both been described separately but using different approaches. Working with Swiss colleagues, I was able to analyse the epidemics in exactly the same way without having to share sensitive data. I found clustering in the UK to be much higher at relaxed thresholds than in Switzerland (34% vs 16%) indicating that the UK database is more likely to capture transmission chains. Down sampling revealed that this pattern is driven by the larger size of the UK epidemic. At tighter cluster thresholds, the epidemics were very similar. I next use these tools to analyse the spread of emerging subtypes A1, C, D and G in the UK. I found both risk group and cluster size to be predictive of cluster growth, which I tested using simulations and a GLM. Growth of MSM and crossover clusters was significantly higher than expected for subtypes A1 and C, indicating that crossover from heterosexuals to MSM has contributed to their expansion within the UK. Numbers were small for subtypes D and G but the proportion of new diagnoses linking to MSM and crossover clusters was similar to A1 and C, suggesting that the same pattern may be emerging for D and G. I conclude by evaluating the accuracy of a method previously described by our group to generate transmission networks from HIV sequences. The interpretation of clustering patterns from phylogenetic trees is difficult because of the absence of a standardised statistical framework. In contrast, a body of work exists that relates disease transmission to networks. Using large simulated datasets, I developed algorithms which eliminate improbable links. I then reconstructed improved UK transmission networks for subtypes A1, B and C and compare network metrics (such as the degree distribution) between risk groups. Together with other evidence, this thesis demonstrates that the UK HIV epidemic continues to be driven by transmission among MSM. The UK epidemic is no longer compartmentalised and the crossing over of subtypes across risk groups has been facilitated by MSM also having sex with women

HIV transmission networks among transgender women in Los Angeles County, CA, USA: a phylogenetic analysis of surveillance data.

BackgroundTransgender women are among the groups at highest risk for HIV infection, with a prevalence of 27·7% in the USA; and despite this known high risk, undiagnosed infection is common in this population. We set out to identify transgender women and their partners in a molecular transmission network to prioritise public health activities.MethodsSince 2006, HIV protease and reverse transcriptase gene (pol) sequences from drug resistance testing have been reported to the Los Angeles County Department of Public Health and linked to demographic data, gender, and HIV transmission risk factor data for each case in the enhanced HIV/AIDS Reporting System. We reconstructed a molecular transmission network by use of HIV-TRAnsmission Cluster Engine (with a pairwise genetic distance threshold of 0·015 substitutions per site) from the earliest pol sequences from 22 398 unique individuals, including 412 (2%) self-identified transgender women. We examined the possible predictors of clustering with multivariate logistic regression. We characterised the genetically linked partners of transgender women and calculated assortativity (the tendency for people to link to other people with the same attributes) for each transmission risk group.Findings8133 (36·3%) of 22 398 individuals clustered in the network across 1722 molecular transmission clusters. Transgender women who indicated a sexual risk factor clustered at the highest frequency in the network, with 147 (43%) of 345 being linked to at least one other person (adjusted odds ratio [aOR] 2·0, p=0·0002). Transgender women were assortative in the network (assortativity 0·06, p<0·001), indicating that they tended to link to other transgender women. Transgender women were more likely than expected to link to other transgender women (OR 4·65, p<0·001) and cisgender men who did not identify as men who have sex with men (MSM; OR 1·53, p<0·001). Transgender women were less likely than expected to link to MSM (OR 0·75, p<0·001), despite the high prevalence of HIV among MSM. Transgender women were distributed across 126 clusters, and cisgender individuals linked to one transgender woman were 9·2 times more likely to link to a second transgender woman than other individuals in the surveillance database. Reconstruction of the transmission network is limited by sample availability, but sequences were available for more than 40% of diagnoses.InterpretationClustering of transgender women and the observed tendency for linkage with cisgender men who did not identify as MSM, shows the potential to use molecular epidemiology both to identify clusters that are likely to include undiagnosed transgender women with HIV and to improve the targeting of public health prevention and treatment services to transgender women.FundingCalifornia HIV and AIDS Research Program and National Institutes of Health-National Institute of Allergy and Infectious Diseases

Lack of Effectiveness of Antiretroviral Therapy in Preventing HIV Infection in Serodiscordant Couples in Uganda: An Observational Study.

BACKGROUND: We examined the real-world effectiveness of ART as an HIV prevention tool among HIV serodiscordant couples in a programmatic setting in a low-income country. METHODS: We enrolled individuals from HIV serodiscordant couples aged ≥18 years of age in Jinja, Uganda from June 2009 - June 2011. In one group of couples the HIV positive partner was receiving ART as they met clinical eligibility criteria (a CD4 cell count ≤250 cells/ μL or WHO Stage III/IV disease). In the second group the infected partner was not yet ART-eligible. We measured HIV incidence by testing the uninfected partner every three months. We conducted genetic linkage studies to determine the source of new infections in seroconverting participants. RESULTS: A total of 586 couples were enrolled of which 249 (42%) of the HIV positive participants were receiving ART at enrollment, and an additional 99 (17%) initiated ART during the study. The median duration of follow-up was 1.5 years. We found 9 new infections among partners of participants who had been receiving ART for at least three months and 8 new infections in partners of participants who had not received ART or received it for less than three months, for incidence rates of 2.09 per 100 person-years (PYRs) and 2.30 per 100 PYRs, respectively. The incidence rate ratio for ART-use was 0.91 (95% confidence interval 0.31-2.70; p=0.999). The hazard ratio for HIV seroconversion associated with ART-use by the positive partner was 1.07 (95% CI 0.41-2.80). A total of 5/7 (71%) of the transmissions on ART and 6/7 (86%) of those not on ART were genetically linked. CONCLUSION: Overall HIV incidence was low in comparison to previous studies of serodiscordant couples. However, ART-use was not associated with a reduced risk of HIV transmission in this study

Transmission of non-B HIV subtypes in the UK is increasingly driven by large non-heterosexual clusters

BACKGROUND: The United Kingdom human immunodeficiency virus (HIV) epidemic was historically dominated by HIV subtype B transmission among men who have sex with men (MSM). Now 50% of diagnoses and prevalent infections are among heterosexual individuals and mainly involve non-B subtypes. Between 2002 and 2010, the prevalence of non-B diagnoses among MSM increased from 5.4% to 17%, and this study focused on the drivers of this change. METHODS: Growth between 2007 and 2009 in transmission clusters among 14 000 subtype A1, C, D, and G sequences from the United Kingdom HIV Drug Resistance Database was analysed by risk group. RESULTS: Of 1148 clusters containing at least 2 sequences in 2007, >75% were pairs and >90% were heterosexual. Most clusters (71.4%) did not grow during the study period. Growth was significantly lower for small clusters and higher for clusters of ≥7 sequences, with the highest growth observed for clusters comprising sequences from MSM and people who inject drugs (PWID). Risk group (P< .0001), cluster size (P< .0001), and subtype (P< .01) were predictive of growth in a generalized linear model. DISCUSSION: Despite the increase in non-B subtypes associated with heterosexual transmission, MSM and PWID are at risk for non-B infections. Crossover of subtype C from heterosexuals to MSM has led to the expansion of this subtype within the United Kingdom

Analysis of the History and Spread of HIV-1 in Uganda using phylodynamics

HIV prevalence has decreased in Uganda since the 1990s, but remains substantial within high-risk groups. Here, we reconstruct the history and spread of HIV subtypes A1 and D in Uganda and explore the transmission dynamics in high-risk populations. We analysed HIV pol sequences from female sex workers in Kampala (n = 42), Lake Victoria fisher-folk (n = 46) and a rural clinical cohort (n = 74), together with publicly available sequences from adjacent regions in Uganda (n = 412) and newly generated sequences from samples taken in Kampala in 1986 (n = 12). Of the sequences from the three Ugandan populations, 60 (37.1 %) were classified as subtype D, 54 (33.3 %) as subtype A1, 31 (19.1 %) as A1/D recombinants, six (3.7 %) as subtype C, one (0.6 %) as subtype G and 10 (6.2 %) as other recombinants. Among the A1/D recombinants we identified a new candidate circulating recombinant form. Phylodynamic and phylogeographic analyses using BEAST indicated that the Ugandan epidemics originated in 1960 (1950-1968) for subtype A1 and 1973 (1970-1977) for D, in rural south-western Uganda with subsequent spread to Kampala. They also showed extensive interconnection with adjacent countries. The sequence analysis shows both epidemics grew exponentially during the 1970s-1980s and decreased from 1992, which agrees with HIV prevalence reports in Uganda. Inclusion of sequences from the 1980s indicated the origin of both epidemics was more recent than expected and substantially narrowed the confidence intervals in comparison to previous estimates. We identified three transmission clusters and ten pairs, none of them including patients from different populations, suggesting active transmission within a structured transmission network

Non-disclosed men who have sex with men in UK HIV transmission networks: phylogenetic analysis of surveillance data.

BACKGROUND: Patients who do not disclose their sexuality, including men who do not disclose same-sex behaviour, are difficult to characterise through traditional epidemiological approaches such as interviews. Using a recently developed method to detect large networks of viral sequences from time-resolved trees, we localised non-disclosed men who have sex with men (MSM) in UK transmission networks, gaining crucial insight into the behaviour of this group. METHODS: For this phylogenetic analysis, we obtained HIV pol sequences from the UK HIV Drug Resistance Database (UKRDB), a central repository for resistance tests done as part of routine clinical care throughout the UK. Sequence data are linked to demographic and clinical data held by the UK Collaborative HIV Cohort study and the national HIV/AIDS reporting system database. Initially, we reconstructed maximum likelihood phylogenies from these sequences, then sequences were selected for time-resolved analysis in BEAST if they were clustered with at least one other sequence at a genetic distance of 4·5% or less with support of at least 90%. We used time-resolved phylogenies to create networks by linking together nodes if sequences shared a common ancestor within the previous 5 years. We identified potential non-disclosed MSM (pnMSM), defined as self-reported heterosexual men who clustered only with men. We measured the network position of pnMSM, including betweenness (a measure of connectedness and importance) and assortativity (the propensity for nodes sharing attributes to link). FINDINGS: 14 405 individuals were in the network, including 8452 MSM, 1743 heterosexual women and 1341 heterosexual men. 249 pnMSM were identified (18·6% of all clustered heterosexual men) in the network. pnMSM were more likely to be black African (p<0·0001), less likely to be infected with subtype B (p=0·006), and were slightly older (p=0·002) than the MSM they clustered with. Mean betweenness centrality was lower for pnMSM than for MSM (1·31, 95% CI 0·48-2·15 in pnMSM vs 2·24, 0·98-3·51 in MSM; p=0·002), indicating that pnMSM were in peripheral positions in MSM clusters. Assortativity by risk group was higher than expected (0·037 vs -0·037, p=0·01) signifying that pnMSM were linked to each other. We found that self-reported heterosexual men were more likely to link MSM and heterosexual women than heterosexual women were to link MSM and heterosexual men (Fisher's exact test p=0·0004; OR 2·24) but the number of such transmission chains was small (only 54 in total vs 32 in women). INTERPRETATION: pnMSM are a subgroup distinct from both MSM and from heterosexual men. They are more likely to choose sexual partners who are also pnMSM and might exhibit lower-risk sexual behaviour than MSM (eg, choosing low-risk partners or consistently using condoms). Heterosexual men are the group most likely to be diagnosed with late-stage disease (ie, low CD4 counts) and non-disclosed MSM might put female partners at higher risk than heterosexual men because non-disclosed MSM have male partners. Hence, pnMSM require specific consideration to ensure they are included in public health interventions. FUNDING: National Institutes of Health

Genomic epidemiology of a densely sampled COVID-19 outbreak in China.

Analysis of genetic sequence data from the SARS-CoV-2 pandemic can provide insights into epidemic origins, worldwide dispersal, and epidemiological history. With few exceptions, genomic epidemiological analysis has focused on geographically distributed data sets with few isolates in any given location. Here, we report an analysis of 20 whole SARS- CoV-2 genomes from a single relatively small and geographically constrained outbreak in Weifang, People's Republic of China. Using Bayesian model-based phylodynamic methods, we estimate a mean basic reproduction number (R 0) of 3.4 (95% highest posterior density interval: 2.1-5.2) in Weifang, and a mean effective reproduction number (Rt) that falls below 1 on 4 February. We further estimate the number of infections through time and compare these estimates to confirmed diagnoses by the Weifang Centers for Disease Control. We find that these estimates are consistent with reported cases and there is unlikely to be a large undiagnosed burden of infection over the period we studied